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Phe} ABSTRACT 


A study of Dynamic Systems with three variable para- 
meters is made by an initial scanning of the basic geo- 
metric properties on the three dimensional Space generated 
by these parameters. From these geometric properties, a 


Root-locus technique that i oe S greatly the amount of 
ecige 


work in analysis and design is developed. This technique 
is extended to systems with k variable parameters. 
Finally, singular lines in the parameter plane are treated 
as a special case of the parameter Space, and formulae are 
derived for a fourth order system leaving the field open 


for systems of higher order. 
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1. INTRODUCTION OF THE THIRD PARAMETER 

This study takes as a starting point Mitrovic's ap- 
proach to the root-finding problem. Although this is 
essentially a two-variable method which will yield only 
two relations out of the three needed for the three-para- 
meter case, the parametric equations obtained will help in 
understanding the nature of the basic geometric properties 
of the parameter space generated by @, 8, and y, the three 
variable parameters. The fact that it is possible to 
obtain only two relations is justified because a point on 
the s-plane is completely defined by only two conditions, 
and therefore any other possible relations will be linear- 
ly related to the two originated by Mitrovic's* root-find- 
ing method. 

Let us now describe briefly Mitrovic's approach, as a 
base for the present study: 

The characteristic equation for a given system can be 


written in a general form as; 


F(S) =a S +a 
n 


or in a closed form as: 


n 
F(S) = ya = = 0 im 2 
k= 


If now S is expressed as: 


S= Cw Be ju /1-0* Le 


where € and Ww are defined in figure 1.1, which is the re- 
presentation of a point in the s-plane that is defined by 


assigning values to C¢ and Ww. 





Co 
n 


Fig. 1.1 - Representation of a point 
on the s-plane 
Equation 1.3, written in another form yields: 
. — re 
S = Wwcos 89 + 5 sin 9) = &e 1.4 
~1 k 

where @ = cos ~(-¢€). Then S~ can be written: 

sK = us eIk? = “(eos k6@ + 3 sin k8), 15 


from which the Chebyshev functions are obvious and are: 


T,, (¢) = Ceaako = GesiGkeom aC) iG 
3 . =i 
_ Sin k@ _ sin(kcos ~€) 1.7 


ag, (6) sin 6 ~ sin(cos-} C) 


LO 


YN 
I 


k wm (-0) + 3 ase u,(-o) | T'was 


ay ml (1) thy se, | 1.9 


Inserting these in 1.2 and requiring that the reals and 
imaginaries go to zero independently, provides the two 
equations: 


n 
) aw ®(-1)*e, (2) = 0 
k=0 


n 
da (-1) hy, (E) = 0 
k=0 


But T(t) = £(u,(¢) -(u,_, (0) 


and upon substitution, the following equations are obtained: 


n 
)(-2) Kaw,” Caer 0 
k=0 


err: 
n 
| k k = 
)(-1) a, Ww U,(6) = 0 
k=0 
where UL. can be obtained from the recurrence formula: 
eee UAE) Ulla) 
eee 
where U,(o) = 0, and U, (¢) el 
Now, the coefficients are defined to be: 
a = db, + Bc, eae tT in J eae lg} 


Va 


where @, B, and y are the independent variable parameters 
of the system under analysis. By substituting this ex- 
pression for ay. in equations 1.11, these equations can be 


rewritten as follows: 


aB, + Bc, jy De hee? 
Ll 
AB. + Bc., + Ye Tae 
where 
n n 
Ee lke k Wa «| k 
By = 2 (-1) "by, TU, (0) By = (1) "by, & “UY (0) 
k=0 k=0 
na Nn 
7 aie k 2 v(- k k 
= Je, ou) Gl eee ce, wm iu (et) 
k=0 k=0 
1s Ji 
Nn Nn 
ee k Paid k 
D, = p(-l)"a, wu, 1 (0) Dd, = )f-1)%a, w vu, (2) 
lan k=0 
Nn n 
ges k eer... k 
Bday /(-1) 2, OU, (6) F, = (-2) "es, @ "0, (0) 
k=0 k=0 


Information concerning the design properties of the system 
must be obtained from these expressions (equations 1.14). 
This information can be; 

a). For a given set of Performance Specifications, 
what values of @, B, and y, the variable parameters of the 
system, are needed in order to obtain the dominance con- 
ditions that will allow us to meet the given performance 


specifications? 


2 


b). The physical meaning of the values of q@, 8, and 
yy, whether or not these values are actually attainable, 
and if not, what modifications must be introduced to the 
system, such that it is possible to obtain the desired 
values. 

c). Existence of dominance conditions, if these 
exist, and the locations of the remaining roots on the 
s-plane, and their influence in the behavior of the system. 

Tf ¢ and ware given, then the set of equations 1.14, 
will contain three unknowns: @, 8, and y. Therefore the 
surfaces generated on the @, 8, y coordinate system can be 
expressed only in a parametric form, generating sets of 
parameter-planes in @ and 8 for a given value of y; or it 
is possible to obtain a three dimensional representation 


which is not a practical tool for engineering purposes. 
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2. REAL ROOTS 
The characteristic equation for a system can be re- 


written as 
OA(S) + BB(S) + yc(S) + D(S) = 0 201 


By making S = - Oo A(S), B(S), C(S) and D(S), become real 


ee 
numbers A, B, C and D, so that the above equation becomes 


OA + BB + 7C + D=0 20 


which is the equation of a plane on the system of coordi- 


nates generated by @ , 8, and y, with intercepts as 


follows: 
on the @ axis @ = - D/A 
6n the P¥axis 8B = - D/B 2.3 
on the y axis y= - D/C 


This is, then the first set of Known surfaces that can be 
used in the three dimensional space generated by a, 8, and 
Y. However,there is need for a graphical representation 
which allows the determination of the desired information 
without tedious use of computarized data, with niles it is 
hard to get a clear vision of the trend of the different 
values for changes in the parameters of the system. In 
order to obtain this graphical representation, it is possi- 
ble to use descriptive geometry as a useful tool for the 
three-dimensional analysis. 

The construction rules are simple and easy to retain 


and with some practice it can be a tool that can be used 
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very often in this study. 

The first step is the representation of a plane ina 
three-dimensional space: A plane is completely defined if 
two lines belonging to such a plane are defined. It is 
possible to draw these two lines by determining the inter- 


cepts of the plane which represents S = - o, in the q@, 8, 


1 
yY space, with the @B and ay planes. (See figure2.la which 
is a representation of the plane S$ = - 0, on the first 
octant of the coordinates axes). Now, by rotating the y 
axis in the plane By in the manner shown in figure 2.la, 
and in amount such that the three axes lie in the plane 
generated by the axes @ and 8, figure 2.1b is obtained, 
which contains all the information concerning the three- 
dimensional graph. However, this time it is possible to 
work with real values and to plot all the information in 
one plane. Of course, once it is learned how to work with 
this construction method it will not be necessary to have 
the three-dimensional graphs. Let us assume we are to 
investigate whatever values of @ and 8, are obtainable 
with given values of § = - a1 and y = Y1° In the three- 
dimensional picture, we would pass a plane parallel to 

the plane @B, through the given value of y = ¥, in order 
to obtain the plane PNM. As can be seen, this is nota 
very practical method even if there exists the possibility 
of building an actual three-dimensional representation of 


the plane S = - O1- By the construction method proposed, 


what has to be done is to project the chosen value of y 


is 








goes OuGbeNGn « 
i ae ; 
-~DA x 
p 
“~ 
reg. 2.15 
Fig. 2.la - Representation of a plane in three dimensions 


Fig. 2.1lb - Representation of a plane uSing projections and 
all coordinate axes in one plane 
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onto the line ab, determine the point N, and project this 
point N onto N', and finally draw a parallel line to ac 
through N'. Then the line de will be the actual represen- 
tation of the line MN of the three-dimensional graph, from 
which we can read the values of @ and 8. This construction 
method can be more easily understood by the use of a simple 
example: 

Example 2.1: Let us assume that a system has the 
following characteristic equation 


st4 (30448)S°+ (604+38)S7+ (1LORE50Y)S + (5Q45y)= 0 2.4 


which can be rearranged in the form of equation 2.2 as 


follows: 
ss 3 2 S 2 4 _ 

G(S) = (3S~+6S°+5)@ + (4S°4+3S°+10S) B + (50S+5)y+S~ = 0 
If S = -5 is to be investigated, then for this value of S: 
G(S) = -220@ - 4758 - 245y + 625 = 0 
ead : 

_ 625i... 
the @ intercept = - =390 = 2.84 

_ 62d 
the 8 intercept = - 475 = 1. SHkS 

# . ae _ 

and the y = 545 = 2.255 


Figure 2.2 is the representation of the positive axes @, §8, 
and y. The negative axes can also be identified by label- 
ing the $8 axis, as the -y axis, the y axis as the -§8 axis. 
Lines ab and ac are the intercepts of the plane S$ = -5 


with the @B and @y planes; dotted lines are extensions of 


I 





a chet ir ge _ 


¥ or(-/3) 


Fig. 2.2 - S = -5 lines in the @By space for Example 2.1 
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the lines on the negative side of the axes. Let us assume 
that we wish to investigate what values of @ and 8B are 
suitable for a given value of y, say y = -5. By setting 
y = -5 and projecting this value of Yy onto line ac, point 
dis obtained from which a line parallel to ab is drawn. 
This is line de, and this line contains the information 


about the values of @ and B that will satisfy equation 2.4 


for S = -5 and y = -5. As a check, let us choose the point 
f; at this point @ = +5, B = +1.58 and y = -5, which sat- 
isfies equation 2.4 if S = -5. 


Another interesting feature of this construction me- 
thod is that due to the linearity properties, once the 
"advance" over the @ axis for a given change in y is deter- 
mined, a grid can be easily drawn without the need of pro- 
jecting each value of y, by just advancing equal distances 
on the @ axis and drawing lines parallel to ab. The loci 
for y = -5.0, y = -7.5, y =+2.5 and y = 5.0 were drawn by 
this method. Similar results can be obtained by determining 
the "advance" in the B axis. In summary, all the necessary 
information is: the intercepts of the plane with the axes, 
the "advance" for a desired increment of one of the para- 
meters. Then the surface generated by S = -5 will be com- 
pletely described. Therefore, for a given value of S = -O, 
the graph will look like the one in figure 2.3. When 
several values of 0 are needed, confusion can be avoided by 
identifying each pair of lines by the common point on the 


@®@ axis (See figure 2.3). 


ES) 





Vv 


Fig. 2.3 - Grid of loci for several values of real roots 


in the aBy space. 
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This method of construction will be useful in cases in 
which it is desired to analyze well damped systems and 
also aS an easy way to construct a grid of values of real 
roots, which can be superimposed over a previously deter- 
mined set of constant € curves, in order to determine the 
possible values of real roots expected for given values of 
oo Boy, GS and WwW. 

Example 2.2: Consider a voltage regulator. 
The compensation schemes have been worked out and this 
example will be analyzed through this study so we can ob- 
tain a better picture of the use of this construction me- 
thod. 

After manipulation of the block diagram the character- 
istic equation of the system is: 

4 6 


S° + (2x10 K 
¢) 


91350) S$” + (3x10°K +4x10°xK 


1 10 
K, +2107 °K 96x01 0°)GW a) abe LO 


#3\e5k20° ys 


TT 
+ (or K K +2.1x1029 = 0 
ap 


Dpll 


iy fomerig K.. = &, Ki, = B, KK = y, and rearranging the 


io ap 


characteristic equation in the form of equation 2 2, it 
becomes: 


nace 8 


(2x10°s?4+3x10 10 


ee 


5*+1019S) a@ + (4x10%s74+2x101%s) B + 


10 4. 3) 5 


2x10) 9 y+974+1350S°+3.85x10 6 


Ee aril ao nea 


$7+36x10 
For S = -1 the equation of the plane is: 


-9702@ - 196008 + 20,000y + 20964.38 = O 


PiM 


20,964.38 _ 


where the @ intercept = - =9702 Zao, 
the B intercept = - ee 0,6 
-19,600 eid 
a alll OAT 
and the y intercept = %,o00' 1.04, 
For S = -5, the equation of the plane is: 


-4275a@ - 9000B + 2000y + 2082.95 = 0 


2032759'5 


the @ intercept = - ~4275 > 0.488, 
the B intercept = - 2084.22 = 0.231, 
and the y intercept = - 20822 — -1.04. 


The results of this problem can be easily checked on 
figure’2.4 forthe Gase lét us take say @ = 0.5," = 0.5 
by drawing a line parallel to ab through the point (@ = 0.5, 


B = 0.5), we get point c on the @ axis, and then by project- 


ing this point onto the line corresponding to o = -1.0 on 
the plane ay, we get y = -0.315 these values satisfy the 
characteristic equation of the system for S = -1.0. Another 


variation can be obtained by selecting a value of y, this 
is done for o = -5.0 in which y = 1.5 was selected, by pro- 
jecting and drawing a parallel line to de we obtain the 
line labeled (o0 = -5.0, y= 1.5). Any pair of values for @ 
and B on this line will satisfy the characteristic equation 
for the system. We have picked the point B, which yields 


® = 2.0 and 8 = -0.385 which checks. 


Z2 


> or (-6) 


\ Sy 
\ ~ “os 


~ 0, 





Fig. 2.4 - S = -5.0 and S = -1.0 loci for Example 2.1 
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3. COMPLEX CONJUGATE ROOTS 

Usually dominance of a certain pair of complex roots 
is desired in the system. By dominance it is understood 
that the system will behave in such a manner that its dy- 
namic performance can be approximated by a second order 
system; that is, that the influence of the remaining 
roots of the characteristic equation will be very small 
on the transient behavior as compared with the influence 
of the pair of dominant poles. This feature will allow 
the designer to obtain a certain degree of control towards 
the achievement of some performance specifications such as 
settling time, natural frequency and damping coefficient 
(overshoot or undershoot). For this analysis equations 


1.11 can be rewritten as; 


k 


2 Ye 1)*bow* u(t) + 2 Y(- 1) Ke wu, 5 (0) 
k=0 
n k 
ty )(-1)¥a,w ku, ,(t) + )(-1)%e,w*u_,(t) = 0 
k=0 k=0 
oral 
si 
a )(-1)*bw ku (t) + 8 ve L) cw *u 0+ 1 )t- ifaw Ku, () 
k=0 k=0 k=0 
a 
+ )(-1)Se,w *u,(£) = 0 
k=0 
where ay. has been replaced by 
a, = Ub, + Bc, + yd. + e,. 32 
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By assuming given values of € and w for a desired domi- 


nance condition, equations 3.1 can be written as: 


n n n n 
a Peed, + B ) GAB act )a4,D,, a )e,,Dy, = 0 
k=O k=0 k=0 }:=0 
n n n n js 
a eee + B Jeae. +y 4,5, am Pe, Ey = Q 
k=0 k=0 k=0 k=0 
_ kk 
where dD. = 5) — U,,_; (0) 
Bia 
P kok 
and E, = (-1) we U,, (C) 


By evaluating the quantities under the summation sign these 


two equations can be further reduced to: 


O, and 


aB, + Bc, Haldia 


ll 
© 


OB. + Bc., + yD, + E, 


These are the equations of two planes in the space aq, 8, y. 
The locus of the values of @, B and y which satisfy the 
desired values of € and Ww a the straight line defin- 
ang the intercept of the two plhanese Et ve known that a 
real root can be represented as a plane in the space a@ 8B y, 
sO now we can determine the intercept of the real-root 
plane with the straight line defined by the intercept of 
the planes described by equations 3.5. There are two 
possibilities: a) the line is not contained on the plane 
and then the intersection will be a point in the space a 8 


Y¥; or b) the line is contained on the plane and then we 
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have a set of values of @ B y on the line that will , 
satisfy the desired conditions. In other words, it will 
be possible to solve a set of three equations with three 
unknowns as follows: 


From equation 3.5 


OB, + Bc, + yD, + E, = 0 
826 
OB, + Bc., + Dee) - 0 
and from equation 2.2 
OA + BB + 7C + D=0 
By Cramer's rule: 
mi -E, (C,C-BD.) + E, (C,C-BD, ) - D(C, D,-C.D,) 
ass K a 
, B, (-E,C+DD,) = Bo (-E, C+D D,) == A(-E,D,+E5D),) 3.7 
—— ES 
and 
. A 


where A= B, (C,C-BD,) = B. (C,C-BD, ) + A(C,D,-C5D)) 


Let us now illustrate this use of Known surfaces in the 
a B y space by an example: 

Assume a third-order system with the following 
charadGterrstic equation: 


= “+ (3044 B) $7 + (40+38)S + 24 + y = 0 Aus: 


Assume now that we wish a pair of complex poles with 


second-order dominance at S = -5 i 8.663 and a third pole 
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at S = -50. What values of @, B and y must we use to 


obtain the desired pole values? 


» 


By applying S = -50.0 to equation 3.8 we obtain: 
7,302@ + 10,0008 - 149y - 125,000 = 0 3.9 
and by substituting §S = -5 Mf 8.66} in equation 3.8 we 
obtain: 


298a¢ + 4008 - vy - 1,000 = O 


and 260@ + 4008 - 30y¥ = O 


By solving the system of equations formed by 3.9 and 3.10 


we get values of 
@ = 7200, B = -5,385 and y = -9,400 


By substituting these values in equation 3.8, we get 


3 2 


S~ + 60S“ + 600 + 500: =0 

or (en Sass os + 100) = 0 Saal 
f) 

or (S + 5)(S + 5 - 8.665) (S + 5 + 8.663) = 0 


We see that we have three degrees of freedom and therefore 
have the opportunity of choosing three of the roots of the 
Pe rae teristic equation for systems of order higher than 
three, once we have made our choice of three of the roots. 
The remaining roots will depend on the values of a, 8B and 

y (except in the very special case of situatim ») previously 
described, in which we can find another root by determining 
the intersection of the line with’ the plane generated by 
this extra root). Therefore, we have the freedom of choos- 


ing three of the roots of the polynomial. However, for 
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polynomials of a degree higher than three, this method does 
not tell us anything about the remaining roots of the poly- 
nomial, and once the choice has been made, the remaining 
roots will be determined. So far, as in the case of the 
parameter space, the designer has the freedom of choosing 
several of the roots of the characteristic equation, but 

1S unable to determine whether or not the remaining roots 
will help in design, or even make the system stable. This 
is a tedious trial and error process, and it is very likely 
that for the chosen roots there will not be a satisfactory 
set of remaining roots. There are some useful methods like 
the Polak method and the Ross-Warren method, employing 
filter compensation, but they are essentially one-variable 


parameter methods and are not applicable in a general sense. 


28 


4. CONSTANT € SURFACES 

Tt was shown in Section,3. thatthe Locus.on the a@ 8 y 
space for given values of © and w was the line defined by 
the intersection of the two planes represented by equations 


3.5 which are: 


I 
CO 


OB, 2. Bc, + yD, + E, 


and OB. + Bc. + YD 


; me = 0 


2 Z 


The above statement can be interpreted also as "for a given 
C¢, the constant w loci on the a B y space are straight 
lines", so that by keeping € constant for each value of we 
there will be a set of two planes in the form of equations 
3.5 with the intercept a straight line in the three dimen- 
Sional space @ B y. By projecting this straight line onto 
one of the planes generated by the coordinate axes we will 
be able to summarize all the information pertaining to this 
straight line in one plane. Figure 4.1 shows the intersec- 
tion of two planes ABC and DEJ as the line GH and its 
extension (dotted line). The projection of GH on the plane 
is the line HF and its extension. As can be seen, equal 
increments in y will produce equal advances along the pro- 
jection. Therefore, linear interpolation can be done 

along the projection line HF, and the information concern- 
ing the line of intersection can be summarized in the @ 8 
plane. The line HF and its extension can easily be deter- 
mined by finding the points H and F. Setting y = 0 in 


equations 3.5, and solving for @ and B we obtain the point 


Pb, 







PROTECTION OF INTERSECTION LINE 
ON Of PLANE 


/ | | 
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Fig. 4.1 - Intersection of two planes in a three 


dimensional space. 
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H, applying Cramer's rule: 


re (-E,C. + E,C,) 
_ AH 
are 
(eee een. ) 
q ae “eal 
and BH = =  y 
where AH = BC. - BoC, 
Similarly for point F, setting B = 0 
a -E,D, + E,D) YS 
AF : 
and cEa—"o 10 


It is possible also, to determine a "measuring-stick”" 


along the line HF by determining the value of ¥Y, given by 
> —A 4.4 


where AF = B.D. = B.D 4.5 


As can be seen, for an increment Ay = Vy we must advance 
along the projection line equal distances HF. We can draw 
the projection line in the @ B plane and by using linear 
interpolation, label along this line the values of y. 
Figure 4.2 illustrates the construction method used in the 
two planes shown in the three-dimensional picture in Fi- 
gure 4.1. For this we have drawn each one of the planes 
represented in equation 3.5 by the method proposed in 
Section 2. The intersection of the two planes in the @ 8 
plane for y = 0 is the point H and the intersection in the 


a B plane through y = Vy is the point F. Therefore, the 
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Fig. 4.2 - Two dimensional representation of the inter- 
section of the two planes shown in Fig. 4.1. 
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projection line is HF. Any other point along this line can 
be determined by: a) picking up the corresponding value of 
eal (see ¥3) and projecting this value onto the @ axis 
through the lines on the @ y plane (ED and BA in this case). 
These projections will determine the points P and Q on the 
@ axis. By drawing a line parallel to BC through P anda 
line parallel to EJ through Q we determine the point 
y= 3 = 3Y\) along the line HF, from which we can obtain 
the corresponding values of @ and 8. b) Linear interpola- 
tion can be used also, by laying off the distance HF along 
the projection of the line for equal increments, Ay = Y1° 
The construction method was used only as a graphical help 
to understand the linear characteristics and the possibili- 
ty of linear interpolation. Everything can be done analy- 
tically. The equation of the projection line can be obtain- 
ed by using the formula for a straight line, given two 
points H and F. An example follows that will illustrate 
Pie wme thod . 

Example 4.1: Assume the system shown in figure 4.3 is 
to have a dominant pair of poles with € = 0.5 and w = 2.0. 


The characteristic equation is 


(Stems) Ge + ys +°5)(S + 10)+ & K_S* + K K,S + K K) 
4.6 
By rearranging terms we get: 
54 + emis + (K eat 73)8° +(K K,+ S225)6S. GK K, + 25)= 0 
4.7 
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(S+t0-5) (S+1)(St5)(S+10) 





Fig. 4.3 - System with three variable parameters, 
Example 4.1. 
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or 


s* - 16.5S°+(a 4 73)S-+(B + S2e SheGe+ ofY ak 25) a= 0 
where @ = K K, B= wk i i 
by setting € = 0.5 and WwW = 2.0 we get the equations: 


44 - y +135 = 0 


4.8 
and 4a -28 + 110 = 0 
by making Y = 0 we get OH = -33.75 
4.9 
and BH = -12.0 
Her P= 0 we get OF = -27.75 4.10 
and y = 24.0 


We can now draw the straight line on the @ 8 plane as 
shown in figure 4.4. 

By inspection of equation 4.7 we can see that one of the 
roots will go in the right-half side of the s-plane if y 
is negative and of absolute value greater than 25; there- 
fore, for the system to be stable y must be y > -25. By 


setting y = 0 the following roots are obtained for equa- 


meron 46 = 
Si2= Pome 1.733 (desired value of roots) 
S, = -14.055 
Sy = -0.446 


The corresponding values of @ and B are -33.75 and -12.0 
respectively. Figure 4.5 shows the location of the roots 


fer several values of Y. Two important features are 
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p STARTING POINT WITH AS FREE paARaAmMETER 
FOR ROOT LOCUS METHOD (SECTION 5S.) 





/ } 
i 


Fag. 4.4 = Stremight line representatiaaggore4, lane 
projection on the @B plane - Example 4.1. 
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noticeable on this graph: a) Due to the relationships 
established in equations 4.8, the desired pair of complex- 
conjugate roots remains in the same place on the s-plane 
for all values of vy. b) The remaining roots of the poly- 
nomial travel through the s-plane in the same fashion as 
on the root-locus. This is essentially a method for deter- 
mining the loci of the roots when there is only one vari- 
able parameter in the system. This fact is justified 
because two degrees of freedom have been used in fixing 
two of the roots of the polynomial. Therefore it is pos- 
sible to obtain a root-locus for the remaining roots of 
the polynomial. The method is outlined in Section 5 and 
proves to be useful in design for systems with several 
variable parameters. Table 4.1 gives the values of the 
remaining roots and the corresponding values of @, B and y 


for the characteristic equation of Example 4.1. 
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Fig. 4.5 - Locus of remaining poles in characteristic 


equation of system in Example 4.1. 
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-~14.055 
BTS 29 
—lglig 7) 

Sg 3 375 


SPST St 225 


-7.25+8.793 


-7.254+10!33 


-7.254+45.33 


Table 4.1 - Corresponding values of S 


system of Example 4.1. 


i), 


-725 2375 
is IS dy 
9 Wome. 795 
27975" 10 . 


-7.25-45.33 


Ba 


S 
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Gand B for 


The straight-line property of the constant Ww lines can 
also be used for interpolation. In order to illustrate 
this let us assume that computer results” yields curves 
Mme = a. y = 7) and Bites = a..y = Yo) on figure 4.6; but 
because of some design considerations the range of admis- 
Sible £iw€ for @ and B is the region OPQR on the @ B 
plane. It is possible then to get a fairly accurate 
picture of the shape of the curves in that region by draw- 
ing straight lines between points on curves A and B having 
corresponding values of Wis and by linear interpolation on 
these lines to obtain curves for intermediate values of y 
between Yy and Yo- Curve C corresponds to y = (y,+ %_)/2 
and was obtained by dividing the straight lines into two 
equal segments. As it can be seen the loop on curve A 


tends to disappear as we approach curve B. 
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Fig. 4.6 - Use of linear interpolation in constant C 


surface using constant we lines. 
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5. MODIFIED ROOT LOCUS 

There are many instances where the performance speci- 
fications for a system require a pair of complex conjugate 
roots for the system. The parameter-space techniques tell 
us what values of the parameters @, B and y are suitable for 
the system to obtain this pair of complex conjugate roots. 
Unfortunately these technigues do not disclose any informa- 
tion about the rest of the roots of the characteristic 
equation for the given values of the parameters unless we 
solve the polynomial, and if dominance or stability condi- 
tions are not satisfied, we will have to go back and make 
another choice of values of @WB y. After several trials it 
1s possible to establish the trend of the roots. However, 
this is a difficult task even if we have a computer program 
available. because in general it is very hard to get a 
clear picture, and in order to get detailed results we 
would need a very complete three-dimensional picture of the 
parameter space. Because of this inconvenience we may use 
a method that gives a picture of the system at one glance 
and tells us something about the remaining roots of the 
polynomial which represents the characteristic equation. 
The name assigned to this method is the "modified root lo- 
cus", because the resulting locus follows the same tech- 
nigues applied to the root locus, except that the desired 
pair of roots remains in the same place on the s-plane 
while the variable parameters @, B and y are changed. 


Let us assume we have a system with characteristic 
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equation G(S, a, B, y), and a pair of complex conjugate 
roots with given values of ¢€ and w is desired. It is 
pesctole) to sceduce thesorder of thescharaeteristic, equa— 


tion by making the division: 


G{S- OF, _y) — a, 8, 5 oa 
S~+ 2Cw + W 
n n 


The guotient of this division Will be a polynomial of 


order n-2: 
Fi(sqwey? By ~y) 5.2 
and the remainder will be 


F,(a, 8, vy)S + Fj(a 8B, ¥) 53 


The division has to be an exact one, in order for 

Se = 20w + w to satisfy the equation G(S, a, 8, vy) = 0. 
Therefore, the remainder (Expression 5.3) has to be equal 
to zero, and for this condition to be satisfied, the values 
of F(a, B, vy) and F(a, B, vy) have to go independently 


to zero. This can be readily seen because for any arbitrary 


value of S equation 5.3 has to go to zero. Then: 


F, (a, may) O 


I 
O 


and F(a, B, y¥) 


This set of equations is, in some cases, the very same set 
of equations represented in equation 3.5,or a linear combi- 
nation of them. Therefore, they represent a line on the 

three-dimensional space @ B y. We can apply the same tech- 


niques used in Section 4 to obtain a graphic or mathematic 
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representation and obtain the corresponding set of values 
of a, B and y that will yield a pair of complex roots 
having values € and Ww. But what about the remaining roots 
of the polynomial? Will they yield the desired dominance 
conditions? Will the system be stable? What are the re- 
quired values of @, B and y? Will these values be phy- 
sically attainable? The answer to these questions lies in 
the quotient of the division represented by Expression 5.2. 
The fact that two of the roots have been fixed means that 
we have implicitly used two degrees of freedom. Therefore, 
there is left only one degree of freedom which is the 
essence of the root-locus method. Its techniques can be 
applied using the following procedure: 

From equations 5.4, we can get any two parameters as 
a function of the third parameter. Let us assume we are to 


solve @ and B as functions of y; then equations 5.4 yield: 
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I 


G, (y) 


and B 


Go(y), 
or, solving for @ and y: 


e4 


G,() 


and Y eae) r) 


or, solving for § and jy: 


so 
| 
a 
el 
= 


and yYy=G (a). 


AA 


Now by substituting the values of say, @ and B in Expres- 


Sion 5.2 (the quotient) we obtain 

moe) 
or if we use @ and y: 

F(S, B) 5.6 
or if we use 8 and yy: 

F(s, @) 


Now we are ready to apply the root-locus techniques to the 
polynomial F(S). Since this polynomial has also to satis- 


fy the characteristic equation we can say that: 
Bis, yy) = 0 Sat 


by rearranging terms of this polynomial we get 


Me MSD re, Sb ogy 5 8 


The root locus for y will be the locus of the remaining 
roots of the characteristic equation. From this root locus 
it 1s possible to establish stability and dominance condi- 
tions. Choose the more convenient value of y, and by sub- 
Stituting it into equation 5.5 get the corresponding values 
of @ and 8; or by a simple look at the straight-line method 
established in Section 4, obtain the physically attainable 
values of @ and B. This modified root locus will also 
allow the designer to select the best choice of the remain- 
ing roots of the polynomial that will help to attain other 


performance specifications, such as bandwidth and gain. 
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The method described above for three parameters can 
be generalized for dynamic systems with k variable para- 
meters. In this case the designer has (k-1l1) choices of 
the roots and there will be one "free parameter" that can 
be varied at will for application of the root-locus tech- 


niques. The characteristic equation will be divided: 


Gil geile Bo wipe cateeevavere iS uae See 


where the denominator of equation 5.9 corresponds to the 
product of the chosen (k-1) roots. 
The quotient of this division will be a polynomial of 


order n -k+1 
Puts sO Ps VY «ese KS 5.10 


and the remainder will be 


F, (0, B,¥,.---k)S* 


ce - FL _5(%,B,y...-k)S + F,_, (Bry...) 
5... dal 
Again, if the division is going to be satisfied with the 


remainder equal to zero a set of k - 1 simultaneous equa- 


tions are obtained. 


F, (a, a yo... eh =00 
F(a, i >. « we 0 

, 5.12 
F,_9(% B, y ... k) = 0 
F (Q, ve. « coda= 0 


AG 


From equations 5.12, it is possible to solve for (k-1) para- 
meters as functions of a Co parameter and by substitution 


14g Expression 5.10 we get: 
F(S, free parameter) = 0 Sy eS} 


jmne application of reot locus techniques to equatuen 35.13 
will produce the locus of the remaining roots of the cha- 
racteristic equation. These techniques are useful in va- 
rious cases, depending on the applications: 

a). for two variable parameters it is possible to set 
one of the roots of the polynomial ata large negative value 
and obtain a root locus for the remaining roots of the po- 
lynomial. 

b). for three parameters the choice of the dominant 
pair of complex conjugate roots can be made, allowing de- 
termination of the remaining roots of the polynomial. It 
is possible to give two of the real roots large negative 
values and determine the remaining roots from the root locus, 
choose the ones that help us to get as close as possible to 
the desired performance specifications. Similar considera- 
tions can be made for k variable parameters, 

The calculations involved are simple and provide an 
overall picture of the system and its possibilities. In 
many cases the quotient is a function of only one of the 
parameters; therefore, the equation for the root locus is 
available without the need of substitutions (see Problem 
5e2). Wath a computer program for root locus the job is 


less laborious. 
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Several examples are available in order to clarify the 
procedure and the techniques. Each one of the cases out- 


lined above is studied. 


Example 5.1: 


+ 7 100 
(S+2Q (S+15) (S+10)(S+ 5)(S+p) 


Kas + KtS 


System with three variable parameters - Example 5.1. 


Fig. So 

The system on Figure 5.1 is to be compensated using tacho- 
meter feedback and accelerometer feedback. We wish to in- 
sure dominance of a pair of complex conjugate roots with 

¢ = 0.5 ald w= 2.0. We wish also to know what are the 
permissible values of the remaining roots of the character- 
istic equations and, if it is possible, to obtain dominance 
conditions the required values of p, k, and k,- The trans- 
fer function for the system on Figure 5.1 is: 


a 100 
57+ (04-50) $4 (500-+830)S°+ (8300-+100B+6250)S7+ (62500100741 500] 5000042 5) 
Be 


FS 


48 


where p= @Q, Kee B, and Be gy 


ByGivueaing by 5- + 2S + 4 we get the quotient: 


5? + (0448) $* + (48@ + 730)S + (730@ + 1008 + 4598) 


SD eleO 


and the remainder 


(4598@+100y - 2008+2884)S + 12080a@ - 4008 - 18292 
a) all 
In order for the remainder to be zero for arbitrary values 
of S, it 1s necessary that 


4598a~ + 100y - 2008 + 2884 0 


oe 2 


II 


and 1208@ - 4008 - 18292 0 


So far, two degrees of freedom have been used; therefore 
we are left with only one of the parameters variable (free 
parameter). The other two are dependent functions of this 
variable parameter (equations 5.12). 

It is easiest to have @ as the variable parameter, 
because most of the coefficients in the quotient, Expres- 
sion 5.11, are functions of @, and also because § is 
easily expressed as a function of @, in the second of 
equations 5.12. From these considerations: 


g — 3020 - 4573 


100 Sees 


By substituting the value of 8 in Expression 5.10 and 
making it equal to zero, we get: 


ae +(a+48) $7 + (48+730)S + (3750@ + 25) = 0 5.14 
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Now, by applying the root locus technique, 


0(S°+48S4+3750)_ _ _ 


5?+485°+730S425 


>peee@(S.+924 = 5603384) (Sit 24 +956. 3385) 4 § 
(G23 647-4120 377) (St200 84742 oT 7GtO.0s43) 


The root locus for the variable parameter @ is shown in 
figure 5.2. For any point on the given root locus it is 
pacers to determine the required value of @ and the cor- 
responding values of B and y are determined from the rela- 
tions expressed in equation 5.12. The graphical method 
may also be used. The root locus in figure 5.2 shows all 
possible values of the remaining three roots of the poly- 
nomial. For low values of @, the real root will be very 
close to the pair of complex poles that are desired as do- 
Minant. For high values of @, the real root will go far 
in the left of the s-plane and the other pair of complex 
roots will be very close to the generating zero of the 
root locus. As can be seen, the system allows in this case 
an improvement in the bandwidth by choosing the value of @ 
that will make the real root lie in the far left of the 
s-plane and selection of the complex roots to insure a 
dominance condition (real part of the complex roots about 
ten times the real part of the dominant roots). Other 
considerations can be made such as what values of @ will 
yield permissible values of K, and K,, Without saturating 


the amplifiers that will provide these gains, etc. 


SO 


+}60 


4j20 


FIXED PAIR’ OF COMPLEX CONS ROOTS 





Fig. S532" — Mo@ufied root locus for remaining roots of 


characteristic equation for gystemeime Filo. Sel. 


yall 


Example 5.822 


A system has the following characteristic equation: 
G(S) = 5? + 16s" + 40s7 4 (a+10)$7 + (B+50)S + (+25) 
It is desired to have a pair of complex conjugate roots 


with C = 0.5 and w = 2.0. By making the division: 


a West 5 Acs en(aee)s- 4 aeesonoes ys Es 


S84 +28 +4 
we obtain the quotient: 


s? + las* + 8s + (0-62) Seulig 
and the remainder 


(-20+B+142)S + (-40+7+273) = 0 5el7Z 


By making the remainder equal to zero we get: 


—-29 + 8 + 142 


and 


0 


-4¢ + ¥Y + 273 


O 


As is easy to see in this case there is no need for substi- 
tution, since the quotient is a function only of @, However, 
the values of B and y will depend on the value chosen for @, 
from the relationships expressed in equations 5.18. 


Rearranging Expression 5.16 for root locus application: 


s? + 1487 + 8S + (@-62) = 0 519 
Or 
3-5 = - 1 5.20 
oo Pie ee = 62 
Gr 
= 1 5.21 


(S°= 1l74e4GS ace Gams falseole 77) =oa 


BZ 


iiewreot lecus for equation 5.21 1s shown in “Ergure See" 
As can be seen it 1S not possible to obtain dominance con- 
ditions. The system will become unstable for very large 
or very small values of &, In this case, with very few 
steps it was possible to obtain a clear picture of the 
system. It is obvious that the root locus must have the 
same shape whatever parameter is chosen. We can check 
this fact by choosing say, y as our variable parameter. 


From eguations 5.18 the expression for @ is: 


By substitution of the expression for @ in the quotient 
(Equation 5.19) we get 


4s° te 568° $9825 £°y + 9273 —F24e = 0 


OG 
— a 1 

3 2 
4(S~ + 14S” + 8S + 6.25 


If K = vy/4, we get: 
eee 
5 2 
(SX + 14S” + 8S + 6.25) 


Or 


1 S223 


K = ~~ _ 
(S + 13.439) (S + 0.28 + jOmoma) 


The root locus for eguation 5.23 is shown in figure 5.4 and 
it is exactly the same as the one on figure 5.3. The part 
which is not present is that corresponding to negative 
values of y. Had we drawn the root locus for negative 


values of y this part would appear. The points where the 


36 


>| -Isiemered eoizy se » 
: = 7*G oT dilexm UL WeasAs 10 SNDOT 3007 pots Tpome— ¢€°s °Srg 


ca v- 9- e- Oe Za ¥ wi- 
| $ n= Lt a 4 a ee dk 5 
5 | 


' / 


cil / 
ort } 


SIOOY XAaIdWOD AO 7 
Yih GIgIS3q / 









/ 
/ 
/ 
at y, 
/ 
/ 
> SiO0a ONINIVAW 
J 7 -3Y4 FHL dO4 SNVOT 100d 
/ 
/ 
es y; ANWId -S 
/ 
/ 
/ 
/ 

/ 

grt 
Vii 


54 


INVId-S 





"tejomeiead saaz Se 4 - ¢€°G oe 1dliexg ut weaskS 20z SNDOT JOON peTFytpom - y°s “Sta 


315, 


generating poles for the root locus of @ appear are the 
points for which the value of yYy make @ = Q. Similarly 
the points where the generating poles for the root locus 
of Y appear correspond to the value of @ that makes y = O. 
We can make a quick check. Let us determine the value of 
@ at point A on figure 5.3. This value yields @ = 68.25, 
and by substituting in equation 5.22 we find that this is 
exactly the value of @ that makes y = O. 

Example 5.3: 

Let us consider again the example 4.1 in which the charac- 
teristic equation is 


oo ~ ews + (0473) 8° + (B +82.5)S + (y+25) = 0 


by uSing the same desired values of © = 0.5 and we = 2.0 
and dividing by S* ora, the quotient is 


s* 4+ 14.55 + (a+440) 


and the remainder is 
(-2@+B8-55.5)S + (-40+y-135) 
Setting the remainder equal to zero, we have 


-4% + y -135 = 0 


ae 5.24 


-2a + B - 55.5 =0 


By comparing equations 5.24 with equations 4.8, in example 
4.1, it is easy to see that these two sets of equations are 
exactly the same; therefore, equations 5.24 represent the 


straight line that defines a constant w on the constant C 
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surfaces. By making a root locus set up with @ as a free 
pabameter we get: 


ad 


— i 
(Ss +14. 5740) 


5.25 
———————E—— 
(E710. a5) (St3-5705) — 


The root locus is shown as the slashed lines in figure 5.5. 
Poles marked D are the generating poles using @® as a free 
parameter. Had y been chosen as the free parameter the 
generating poles would be the ones marked as B on figure 
5.5, and the root locus using Y as the free parameter would 
be the part shown as heavy lines, plus the slashed lines. 
For roots at points between B and D the value of @ would be 
negative. This fact can be checked in figure 4.4 where the 
corresponding points B and D are marked. The point of im- 
portance in this case is that since root locus techniques 
consider only positive values of the variable parameter, 
for each parameter chosen as free parameters, we will get 
only part of the overall root locus. Therefore, if the 
results with a positive value of the free parameter are not 
satisfactory, it is good practice to make a quick check for 
negative values of the free parameter. This can be done 
easily by starting at the generating poles and going in the 
opposite direction. For this specific example, we can see 
that for negative values of @, the root locus will go from 
D towards Bisee Figure 5.3), and eventually one of the roots 


will cross the jW axis and the system will become unstable. 
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Example 5.4- Example with parameters product. 

Gaeke ware Some cases where™= theveharacteristic equatron=con-— 
tagns products of parameters. This techhigué is aipplicebite 
Only if the quotient of the division is as,funetionseot only 
one of the parameters. Let us assume a system with charac- 
teristic equation: 


s* + (a415)S° + (200+200)S* + (7+600)S + aB + 6000 


520 
and the system is required to have a pair of complex con- 
jugate roots with € = 0.5 and w = 10. From the division 
we get quotient 


s* + (045)S + (100+50) 5.27 


and the remainder 
(-2000+y7-400)S + (@B-10000+1000) = 0 5.28 


The root locus technique can be applied to the Expression 


57, thus; 


a@(S+10) . 1 
5 —_— 
(S~+5S+50) 


Of 


a(S+10 
(S+2.5-6.6444) (S+2.5+6. 0444) 


and the root locus can be seen in figure 5.6 
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DESIRED PAIL OF COMPLEX CONJUGATE POLES 
(BY COINCIDENCE ON THE ROOT LOCUS FOR THE REMAINING ZOOTS) 
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Fig. 5.6 ~ Modified root locus for Example 5.5 - @ as free 


parameter . 
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The choice for dominance can be made as a pair of poles at 


S = -20. From the root locus we get 
&@ = 35 


and from Expression 5.28: 


B = 1000q@ - 1000 





a'4 
_ 34000 _ 
prox 55 = OF) 
and 
y = 200M + 400 = 7400 


we could also obtain the values by multiplication. 


G(S) + ($420) * ($7+10S+100) = 0 


G(S) = S2coc” aE 90087 + 8000S+40.000 = 0 


By comparison of the two expressions for G(S) (Eqs. 


and 5.29) we get 


o2+ 15 = 50 


Or © @& = 35 
20@ + 200 = 900 
ana yaweoo"= sooo -y ="74oc 
aa dee 60g = 40.000 : 8 = 22822871 


which checks with the results above. 


This is a special case where the quotient (Eq. 


5.27) 


is a function of only one of the parameters. In other 


cases the root locus technique is not applicable because 


the free parameter appears in the quotient as K, K , K, 


ol 


3 


etc., depending on the complexity of the parameter products 
in the characteristic equation. However, there exists the 
possibility that the designer has enough freedom to modify 
the system configuration and make the selection of the va- 
Yriable parameters in such a way that no products appear, 

or in case of products, that the quotient be a function of 
only one parameter. All this can be achieved by changing 


some feedback paths or series compensation schemes. 
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¢. SINGULAR LINES AND SELF-ADAPTIVE SYSTEMS 

In the analysis of systems with two variable parameters 
by application of the transformations proposed by Mitrovic 
(coefficient plane) and Siljak (parameter plane), the pre- 
sence of systems with what is known as singular lines was 
detected. A.summary of the procedure that led to the 
detection of this type of system follows. This leads toa 
better understanding of what is meant by singular lines. 
From Section 1, the characteristic equation for a given 


system can be written as; 
Bes) <= a_S® +ea 

n 
69h 


Repeating the manipulations of the above equation outlined 


in Section 1, the following expressions are obtained 


al 
yr) Kaya* eso aee 
k=0 


©..2 
n 
k k os 
)(-1) aoe UC 
k=0 
Now by defining the coefficients as 
a, = Ob, + Be, + d, e.3 


and the substitution of ay. in Equation 6.2, this set of 


equations can be rewritten as the system of tinear equa-— 


tions 


Il 
© 


aB. + Bc, + D 


1 i} 


aB, + Bc. re Doe 0 


63 


where Bie Bos Cys Coe Dy and D, are defined as in Equa- 
tions 1: l5a When ~vwelues.of CC and ware assigned, and the 
values of B 


Bos Cie Coo D, and Do are calculated, each 


1’ 1 
of the equations 6.4 will represent a straight line on the 
aB plane and in a geometric sense, the intersection of 

these two straight lines will yield the values of @ and 8B 
that satisfy the system of equations 6.4, and insure that 
the characteristic equation will have a pair of complex con- 
jugate roots with the assigned values of ¢ and Ww. But the 
intersection of two lines in a plane is not always a point. 
Two other cases are possible: a) the two lines are parallel 
in the @B plane and in this case no pair of @, B values will 
satisfy the system of equations 6.4 (see Fig. 6.1b), or 

b) the two lines are the same in which case they intercept 
on an infinite number of points and any point along this 
line will yield a pair of complex conjugate roots with va- 
lues € and Ww. The three different cases are shown in Fi- 
gure 6.1. Figure 6.lc shows the case of coincidence of the 
two lines; this is known as a singular line case and its 
name is derived from the matrix algebra manipulation® On 

the system of equations 6.4. The systems with singular line 
characteristics find application in self-adaptive systems 
where changes in the dynamic behavior depends on the ability 
of the system to cope with variations of the parameters. 

The identification of the values of the parameters and its 
changes in these self-adaptive systems is done by observa- 


tion of the output of the system. Identification is a pro- 
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Fig. 6.1b 


o> 


Pee, Onlc 


Different cases of mapping of a point on the s-plane 


in the @B plane. 
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blem beyond the scope of this work, but literature on the 
subject is available (see reference 5 and 6). Once the 
change on the parameters is identified the self-adaptive 
systems use sensitivity models which are fairly good for 
very small changes in the parameters. This is because of 
the approximations made in the selection of the sensitivity 
models. The idea of using systems with singular lines ina 
self-sdaptive system came up because it yields a precise 
mathematical expression (equation of a straight line) that 
relates the changing parameters. This expression ean be 
simulated by analog or digital means, so the problem of 
steadyness of the dynamic behavior of the system is solved 
in part. A sensitivity study on the remaining roots of 
the polynomial will tell if the desired conditions are 
going to be jeopardized by the given changes in the para- 
meters and what tolerances can be allowed in the change of 
the parameters. The attractiveness of the use of systems 
with singular lines in self-adaptive systems lies in the 
fact that it is only necessary that one of the parameters 
be adjusted for changes detected in the other parameter 
and vice-versa. As was said before, the expression regula- 
ting these adjustments is the equation of a straight line, 
and if sensitivity studies are feasible, the tolerance in 
the changes of the parameters can be determined. 

So far, the detection of systems with singular line 
characteristics is,a,difficult,job,, and even inwthe case 


where singular lines are obtained the values of € and we 
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which generate these singular lines might not be satisfac- 
tory for the desired performance specifications. The sub- 
ject is relatively new and research is going on on the 
matter. Therefore, the possibility of an easy method of 
detection and set up of singular lines might appear with 
time. There are computer programs available for detection 
of singular eee es 

The results obtained in Section 4 and 5 of this work 
led to the thought that singular lines on systems with two 
variable parameters can be considered as special cases of 
a system with three variable parameters in which one of the 
parameters is independent of the other two and has a fixed 
value. This value is obtained from one of the planes re- 
presented by equations 3.5, It 1s clear that the above men- 
tioned plane is parallel to the one generated by the axes 
represented by the other two parameters. 

Figures 6.2 and 6.3 will help to clarify the concept 
above stated. Figure 6.2 shows the general case where € 
and Ww have been given and the point corresponding to the 
given value of S maps in the three-dimensional space as 
the line AB which is the locus of the values of @, B and y 
that will produce the desired pair of roots with parameters 
for S : € and won the s-plane. Line AC is the projection 
of line AB onto the @f plane. Figure 6.3 shows the special 
case where plane A is parallel to the @B plane through y = ey 
The intercept of planes A and B is the line AB, that esta- 


blishes the relation between @ and $ as a straight line on 
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Fig. 6.2 - General case of mapping of a point on the 


s-plane into a straight line in the three 


dimensional space aBy. 
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Fig. 6.3 - Special case of mapping of a point on the 


s-plane with generation of singular line. 


ow 


the @ B plane that passes through y = y (Plane A in this 


1° 
case). So much for the geometric interpretation; the ma- 
thematical interpretation in terms of the division method 
outlined in Section 5 is that from one of the two equations 
obtained from the remainder, F, (a, B, y) = Ono F,(@, B, y)= O 
a specific value of one of the parameters is obtained and 
by substitution of this value in the remaining equation, we 
will get the equation of the straight line relating to the 
other two parameters. An example follows. 

Example 6.1: 

Let us assume we have a system with characteristic equation 
G(S) = 57+ (10a+10) $2+ (40045 B+107+30) $°+(800410 B+252) 


+(70@+10 8+25y+340) 6.5 


and also assume that somehow we have found that a pair of 
complex conjugate roots with € = 0.5 and w = 2.0 will 
yield equation F, (a, B, y) = 0 from the set of equations 5.4, 
defining a plane parallel to the @ B plane. Division of 
equation 6.5 by 5? + 25 + 4 yields as quotient: 

5? + (10@ + 8)S + (20@ + 5B + 10y + 10) 6.6 
and the remainder is 

(-20y + 200)S + (-10@ - 10B - 15y + 300) Gay 


By letting the remainder equal to zero we get 


y-10=0 6.8a 
10m + 108 + 15y - 300 = 0 6.8b 


The graphical representation of equations 6.8 is shown in 
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figure 6.4. Now, by continuing, from equation 6.8a: 

y = 10 6.9a 
and from equation 6.8b, substituting this value of y: 

a> B- 15 =0 6.9b 
Now substituting back in the characteristic equation the 
value of y we get 
G(s) = or (10@+10) SP + (40a+5 B+130) 34(800+10B+ 252)S+(70@ +10B8+590) 

€.10 

and dividing by 5? + 2S + 4 we get the quotient 


5? + (100+8)S + (20@ + 58 +110) 6. Bik 


Expression 6.11 is the same as 6.6, had the value of y been 
substituted back in this expression directly. The remainder 


for the above division is: 


-10@ - 10 +»50 


and by making the remainder equal to zero: 
@a+fBP-415=0 6.12 


This is a singular line case in the @ B plane, because by 
assigning to y the value of 10 there are only two degrees 
of freedom left: @ and B in equation 6.10. According to 
what was stated in Section 5 and taking G, (S) of equation 
6.10 as the characteristic equation of the system with two 
variable parameters and using the two degrees of freedom in 
setting the value of a pair of complex conjugate roots 

s = -1.0 * 1.7323, the remainder of the division should be 
a set of two linear simultaneous equations with a unique 


solution for @ and 8. 
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Fig. 6.4 - Singular line case for system represented by 


Equation 6.5 - Example 6.1. 
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Equation 6.12 is the only requirement for the remainder to 
be zero and therefore an infinite number of sets of values 
of @ and 8 can be assigned. The root locus for the remain- 
ing roots of the polynomial is shown in figure 6.5, with @ 


as the free parameter. This was obtained by substitution of 
B=15 -@ (Equation 6.9b) 


reo equation 6. ll; “then: 


s* + (low + 8)S + (15a + 185) =0 6e13 


or 


10a¢(S + 1.5) 


aa ae le 


== 


If 10@ = K then 


K(S4+1.5 
(S+4-313) (S+4+ 313) 


By inspection of the root locus, and with the criterion of 
existence of dominance being that the real part of all of 
the remaining roots of the polynomial must be greater or 
egual to 10 times the real part of se + 2S"+ 4 at canebe 
seen that the root locus enters the desired region at 
point A and a section of it leaves the region at point B. 
Therefore, for values of @ between these two regions we 
have acceptable dominance conditions. For point A the 


value of @is 2.41. By setting @: 
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Fig. 6.5 - Modified root locus @ as free parameter. 


Example 6.1 - Singular line case. 
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we can have tolerances for & of + 33% from the chosen 


value. Then for 8B we get: 


B = 13.19 


with tolerance for B of i 4.65% from the assigned value. 


It has been proven that the singular line in the a 8 
plane can be taken as a special case of the three-parame- 
ter case. Unfortunately, the detection of planes parallel 
to either of the planes containing the coordinate axes is 
a job as difficult as the detection of singular line in the 
a B case. From the analytic point of view it will be pos- 
sible to modify the coefficients of the different powers of 
S in such a way that a plane originally at an angle with 
the planes containing the coordinate axes becomes parallel 
to one of them. But from the practical point of view it is 
difficult to get this modification into the system, because 
these coefficients come from manipulation on the block dia- 
gram and usually affect several terms of the characteristic 
equation. Another limitation may be that the desired values 
are not physically attainable, due to the intrinsic compo- 
Sition of the elements of the system. From preceeding sec- 
tions it has been known that a point on the s-plane maps 
in a straight line in the a 8 y space and this straight line 
has a projection on the @ B plane. Therefore, it might be 
easier from the designer's point of view, to"simulate" a 
Singular line case by using three variable parameters, and 
instead of one set-up for adjustment of parameters, use two; 


that is: for each change in say, Y, it would be necessary 


yee, 


to modify accordingly @ and B. The expressions that will 
determine the modification procedure are precise mathemati- 
cal expressions represented by equations 5.4 and they can 
be simulated by analog or digital means available. How- 
ever, there are considerations such as weight, space and 
economy that will put a burden on a design of this type, 
but in the event that they can be overlooked, we will gain 
the ability of choosing the pair of dominant poles that 
will help in the achievement of the performance specifica- 
tions and obtaining solid tolerance limits on the parame- 
ters. Engineering judgement must be used in this case to 
weight all these considerations. 

When the design considerations stated above become a 
real problem and a singular line case is desirable, it is 
worth while to study the coefficients of the polynomial 
which represents the characteristic equation, and by alge- 
braic manipulation, obtain the singular line case. The 
method outlined in the following section for a fourth order 
polynomial can be used in higher-order polynomials. The 
results obtained are quite satisfactory and a clear picture 


of what has to be done to obtain singular lines is presented. 
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6.1 GENERATION OF SINGULAR LINES IN A FOURTH-ORDER 
POLYNOMIAL WITH THREE VARIABLE PARAMETERS. 
A Fourth-Order characteristic equation of a system 


with three variable parameters can be written in general 


form as: 
oe 4. (b,Q+c. Bte,yt+Ff \s? + (b.Bt+c., Bte,yt+Ff 7 + 
3 3 3 3 2 2 2 2 
6.14 
(b, G+c, Bre, y+£)S + (b atc Bre y+£.) = @ 
By making the substitution of 
q.. = e,Y + fy 65 
equation 6.14 can be written as: 
—_— (b_ a+c. B+d a + (b.a+c. B+d ii + (b. @+c, B+d.)S + 
3 5 3 Z 2 Z 1 al il 
6 Jali 


(b a+c B+d .) 26) 


It is clear that the case of singular lines will depend on 
the configuration of the polynomial, that is, on the coeffi- 
cients of the different powers of S. Therefore, it 1S pos- 
sible to determine what coefficients must be affected in 
order to obtain singular lines. If it is assumed that the 
polynomial has singular lines for values of x = 20w and 

y = a, division of the polynomial by: 5? + x S + y yields 
the follwoing results: 


Quotient: 


2 2 
S~ + (b,a+c.B+d.-x)S 4 (b5-xb,) @ +(c,-xc,) B +(d,-y-xd,+x ) 


ay ll 7 
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The remainder is: 


F,(a,B)S + F, (a, B) 6.18 
where: 
F, (a, 8) = (b.-yb,~xb spas )@ + (c,-yc,-xc ae )B + 
lias eos oe 3s 1? eae” Se 
2 3 6.19 
(d, yd,t+2xy-xd5+x d.-x ) 
and 
F. (a, B) = (b -yb,+xyb,) @& + (c -yc,t+xyc,)B + 
620 


2 Z 
(do -ydoty txyd.,-yx ) 


Singular lines can be generated for each of the following 


Cases: 


O for arbitrary values of @ and 8. Then: 


Case I: F, (a, B) 


b, - yb, - xb, + xb, = 0 \ 6.21 

Gue—Sycy -— xen it x7c = 0 6822 
ala 2 call : 

d. —- yd. + 2x%y - xd, + xd - x? = 0 6.23 
il. 3 2 3 


which is a system with three equations and two unknowns (x 
and y). From inspection of the system the following conclu- 
sions are drawn: 

a). It is possible to solve for x and y by utilizing 
only equations 6.21 and 6.22. 
1! dq and 


d., which in turn are functions of y (equation 6.15). 


b). Equation 6.23 is a function only of d 


Therefore, once the values of x and y are found, Y can be 


set at the value that will satisfy equation 6.23 for the 


7S 


values obtained for x and y. It is worth noting that the 
modifications in equation 6.23 do not affect the values of 


x and y obtained before. Then, from equations 6.21 and 


6.22: 
by - bic 
20w ee Cc b. - boc 
6.24 
2 
2 ’ _ bi - xb, + xX b, 
n Y b, 


and the third equation to be satisfied once x and y have 
been found is: 


+ 2xy - xd, + “7d - “7 = 0 Ome) 


d, - yd 5 3 


1 3 
By substituting dh. = Syy + fy and using the values obtained 
for x and y we can solve for Y. 

If the values of € and ware not satisfactory to the 
designer, it is possible to know exactly what coefficients 
can be modified in order to get the desired values. In this 
case if for example y = w * is satisfactory but we wish to 


Cc b, and C. can be modified at 


ky Qi 3 


will without disturbing the value of y. From the above 


change x = 20w, Enen Cc 


coefficients the designer might be able to choose the most 
convenient methods for the acutal modification, such as 
feedback paths, series compensation devices, etc. The 
above considerations are applicable to cases II and III 


which follow. 


fo, 


Case Il. F., (a, B) = 0 for arbitrary values of @ and B. Then: 


bo - yb. + xyb, = 0 


6 26 


II 
© 


Co 7 Yen + xyc, 


d - yd, + = ate xyd, - wR = 0 


From the first two equations we get: 
boc5 - boc, 


n bo°3 - bc, 


and 


Cl yd. ar a4 + xyd., - wc = 0 6.28 


Again, by substitution of d.. = Qe, + £, and using the 


values obtained for x and y, it is possible to solve for y. 


GasenDIi. F, (a, B) = F. (a, B) for arbitrary values of @ 
and B. Then: 
2 “ 
by ~ yb, - xb, + X b, = bo - yb. + xyb, 
oe ve we Ke, TL re =c - yc, + xyc 6.29 
1 8; 2 3 O 2 2 i 
Gee MG 2 - xd. + 24 - Re = d -yd.~t+ 24. xyd.- eo 
eT es a an) ee ae oe oie YOR37¥ 


From the first two equations we get 


. (b. ~ b.) (cy = co) + (b, = bo) (c., — Co) 
n b,(c, > = cy) + c, (by - by + bo) 
2 6.30 
2 : en - xb, + xX b, a) 
fi =e “eee Boe = 


3 3 2 
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and 


wv) 3 2 2 
3 + 2xy - xd. SEY a i eae d- yd. — y + ayeaj- yx 


de leva 3 


ai 


That with the given values of x and y will yield the value 
er y. 

The important feature of the results stated above is 
that the designer knows exactly what coefficients of the 
characteristic equation affect directly the values of € 
and a and hence, what parts of the actual system must be 
modified in order to obtain the desired values. 

The solutions proposed above are for the case where 
the "free parameter" is y; that is: one of the planes 
considered in Section 6 is parallel to the GB plane through 
the calculated value of y. The quotient (equation 6.17) 
will be then a function of @ and B. From the remainder, & 
can be expressed as a function of 8, or B as function of @, 
and root locus techniques applied in order to determine the 
location of the remaining roots of the polynomial. A sensi- 
tivity study can be made on the root locus. For applica- 
tions in self-adaptive systems there will be need of only 
one set-up for modifications when changes in either of the 
parameters occur. The following example will illustrate 
the technique. 


Example 6.2: 


A system has the following characteristic equation: 


54 + (100+10) S74 (400+5 B+130) $~+ (800+10B+10y+5)S + 


(700+10B+y+200) 6.32 
81 


Question: What values of ¢€ and w are attainable in sin- 
gular line cases and what are the required values of y? 
For Case I, substituting the values of the coefficients 


in equations 6.24: 


y _ ay = 4 or eh) = 2 6.33 
nN al 
Y = 24.7 


Substituting the value of Y in equation 6.32 and dividing 
by 5? + 258 + 4 -we get: 
Quotient: 


a + (100+8)S + (200+58+110) = 0 6.34 


Remainder : 
-10o@ -108 - 215.3 = 0 6.35 


Now root locus techniques can be applied to equation 6.34 
after substituting B as a function of @ from equation 6.35. 
The expressions for the quotient and the remainder can also 
be obtained from equations 6.17 and 6.20 simply by substi- 
tuting the values of the coefficients. 

For Case II, substituting the values of the coefficients 


in equations 6.27: 


er U2h) mel eecteens Sige 1765 

y= w*=0.2 or w = 1.4142 &. 36 
nN nN 

y = 46.5 


Substituting the value of the y in equation 6.32 and divid- 
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iene by 5? twOuwbdS +a2 evemgein: 
Quotient: 
5? + (1004+9.5)S + (3504+584+123.25) = 0 on 37 
Remainder: 


(42.5@¢ + 7.58 + 389.375)S = 0 


For Case III, substituting the values of the coefficients 


in equations 6.30: 


20 wu se =4 '6r Ce 0r 5 


eo? = a 
= i Or we = 1 Susy’ 
y = 204/9 = 22.667 


Division yields 
Quotient: 
5? + (100+9)S + (300+58+120) = 0 eo 39 


Remainder: 


(400+58+102.667)S + (400+58+102.667) = 0 6.40 
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7. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK 

In this work it is shown that the increase in degrees 
of freedom in a system really helps the designer in the 
task of obtaining from it the desired characteristics. 
Mitrovic's approach has proved to be useful in the parame- 
ter plane or when products of parameters are present in 
the characteristic equation of a system. However, when the 
number of variable parameters is three or more it is diffi- 
cult to obtain graphical information. The tabulated data 
obtainable from a computer comes only in quantized form and 
it is difficult for the designer to use it in seeking his 
objective, which is to make the best choice of the variable 
parameters. In the field of singular lines a method of ge- 
neration of singular lines for a fourth-order polynomial 
was established and it is expected that formulae can be 
obtained for systems with characteristic equation of higher 
order. With this method the values of € and w at which 
singular lines are expected and the value at which a third 
parameter must be set are obtained. In addition, the mo- 
difications in the system for which the values of ¢ and Ww 


are changed are indicated. 
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